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We present the results of an analytical theory and simulations of the field inside a cavity created 
in a dipolar liquid placed in a uniform external electric field. The analytical theory shows that 
the limit of continuum electrostatics is reached through a singularity in the microscopic response 
function responsible for a non-decaying longitudinal polarization wave. Fields in microscopic cav- 
ities are much different from macroscopic predictions, and low-polarity dielectrics are predicted to 
have a continuum limit distinct from the solution of Maxwell's equations. Computer Monte Carlo 
simulations never reach the standard continuum limit and instead converge to the new continuum 
solution with increasing cavity size. 
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Electric fields within cavities in uniformly polarized di- 
electrics are commonly calculated by the rules of macro- 
scopic electrostatics In case of an empty spherical 
cavity carved from a polarized dielectric, the connection 
between the field inside the cavity, Ec, and the macro- 
scopic (Maxwell) field in the dielectric, E, is particularly 
simple: 
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Here, e is the dielectric constant and Eg is the uniform 
external field. Equation ^ is widely used forproblems 
related to inserting nonpolar and polar [J impuri- 
ties into dielectrics and, more recently, for electrical and 
optical properties of nanopartices in polar matrices Q. 
More generally, the calculation of fields within cavities 
in systems with dipolar interactions is fundamental for 
the mean-field formulation of susceptibility (magnetic or 
dielectric) of such media [5]. Despite its importance, the 
limits of the applicability of Eq. ^ have never been stud- 
ied. In particular, one wonders at which cavity size the 
laws of macroscopic electrostatics cease to apply and one 
needs to deal with microscopic electric fields. The so- 
lution of Maxwell's equations is sensitive to the macro- 
scopic shape of the dielectric samples Q and, given the 
short length of correlation decay in liquids [y], this pic- 
ture may break down at a microscopic length-scale. 

Equation ^ lends itself directly to tests since it pre- 
dicts two physically significant results: (i) the cavity field 
is independent of the cavity radius, i.e. the cavity length- 
scale does not enter the final result, and (ii) the dielec- 
tric is infinitely polarizable, i.e. the internal field of po- 
larized dipoles within the dielectric screens the external 
field and essentially no field is expected inside a cavity in 
dielectrics with high e. This Letter analyzes these predic- 
tions by using analytical formulation for the microscopic 
dielectric response and Monte Carlo (MC) simulations of 
cavity fields created inside the model fluid of dipolar hard 
spheres (DHS). 

A cavity within a dielectric can be described by ex- 
cluding the polarization field from its volume 0, Q . The 



generating functional of the Gaussian polarization field 
P can then be written as 
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The product of J-functions in Eq. ^ excludes the po- 
larization field from the cavity volume f2o- The asterisks 
between vectors denote both the space integration and 
tensor contraction, and the bath Hamiltonian Hb [P] de- 
scribes Gaussian fiuctuations of the isotropic polar liquid 
characterized by the response function XsO^)- In dipo- 
lar liquids with axial symmetry, this response function 
expands into longitudinal (L) and transverse (T) projec- 
tions [9] 
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where 3^ = kk and J-^ 1 — kk are two orthogo- 
nal dyads and S^'^(k) are the structure factors which 
depend only on the magnitude of the wave-vector fc; 
y — (47r/9)/3pm^ is the standard density of dipoles m 
usually appearing in dielectric theories, p is the number 
density, and /3 is the inverse temperature. 

The constraint imposed on the polarization field to 
vanish from the cavity breaks the isotropic symmetry 
of the system and produces a non-local response func- 
tion x(ki , k2) (2-rank tensor) depending on two wavevec- 
tors Jii]. This function is obtained by taking the second 
derivative of ln[G(A)] in the auxiliary field A in Eq. ([2]) 
and setting A = 0. The result is [lol |: 
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where the correction term x'^°''''(ki, k2) accounts for the 
effect of the cavity excluding the polarization field from 
its volume. 

We now consider a spherical cavity of radius i?o in- 
side a microscopic dielectric liquid and use the response 
function from Eqs. ([3]) and ([4]) to determine the cavity 
field. For a dielectric in the uniform external field Eq, the 
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FIG. 1: Cavity field calculated from Eq. (0 with two cav- 
ity sizes indicated by the distance of closest approach ri — 
Ro/<y + 0.5 in the plot. The points where obtained by numer- 
ical integration in Eq. ((7} with S'"''^{k) from MC simulations 
(ri — 1.0), while the dashed lines refer to the use of analytical 
S^'^{k) from Ref . [Til. The integral is calculated numerically 
before the appearance of the singularity on the real axis [Eq. 
(|11|) ] and by summation over the poles when the singular- 
ity falls on the axis. The two methods give identical results 
when numerical integration is justified. The upper and lower 
solid lines refer to two continuum limits, Eq. (|10|) and Eq. 
(Ul), respectively. The dash-dotted line refers to the lattice 
summation [Eq. (|12|l ] instead of continuous integration in Eq. 
^ taken for a cubic ceU of iV = 108, ri = 2.0. 



projection of the field inside the cavity on the direction 
§0 = Eiq/Eq becomes: 



Ec = Eo+eo-T*x*'Eo-eo. 



(5) 



Here, Eq — EoSkfi is the Fourier transform of the exter- 
nal field and T is the k-space dipole-dipole interaction 
tensor excluding the hard cavity core with the radius of 
closest approach Ri — Rq + <j/2 {a is the hard-sphere 
diameter) : 
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In Eq. jSI), Dk = 3kk - 1, k = k/fc, and jn{x) is the 
spherical Bessel function of order n. After some algebra, 
one arrives at the following equation 
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Equation ([7]) is the central result of our analytical 
model. The first term in Eq. ([7]) is the local Lorentz 
field [2] which appears in our formalism as the field in- 
side small cavities [l^] of the size much smaller than the 
length of dipolar correlations in the liquid. The oppo- 
site limit of macroscopically large cavities turns out to 
be harder to derive. 
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FIG. 2: Shaded area indicates the region of existence of a 
real-fc singularity k* in the plane of the distance of closest 
approach Ri = Rq + a /2 and the solvent dielectric constant 
e. The analytical form of S'^(fe) (Ref. [ij) was used to solve 
Eq. ini at p* = 0.8. 



The macroscopic (continuum) limit corresponds to the 
neglect of the fc-dependence in the correlation functions 
representing dipolar fluctuations of the polar solvent. If 
all the functions S^''^{k) and A{k) in the parenthesis 
under the /c-integral in Eq. ^ are replaced by their cor- 
responding k = values, Eq. ([7]) transforms to Eq. ([T]). 
However, the motivation for replacing A{k) with A{0) is 
not clear since this function [Eq. ([8])] decays on approx- 
imately the same length-scale as j1{kRi) in Eq. ([7|). It 
turns out that, if one employs the identity 
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and applies the "continuum" limit to the term in the 
parentheses, one gets an alternative expression for the 
cavity field 
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The direct /c-integration in Eq. ([7]) shows that the 
actual solution branches between two continuum limits 
(Fig. [T]) through a singularity point which appears when 
a pole of the longitudinal function 



S^{k*) + 2A{k*)=0 



(11) 



falls on the real axis (Im(fc*) = 0). Equation (|10p accu- 
rately describes the cavity field at low polarities switching 
to a solution close to Eq. ([T]) through a singularity. The 
appearance of a real-Zc singularity prevents us from us- 
ing numerical integration. The high-e parts of the plots 
in Fig. [T] have been calculated by residue calculus using 
the analytical representation for S^''T[k) from the mean- 
spherical approximation [13] re-parametrized to give the 
exact fc = limits in terms of e and y [11]. The values of 
e{y) have been taken from MC simulations. 

The real-axis singularity in Eq. ^ signals the appear- 
ance of a non-decaying polarization wave induced by the 
cavity and radially propagating from it through the en- 
tire liquid. This longitudinal polarization wave is ter- 
minated at the boundary of a dielectric sample where it 



3 



3 - 




20 40 60 80 100 



FIG. 3: Cavity field obtained from MC simulations of the 
DHS fluid with a spherical cavity at the center of the simula- 
tion box. The points in the plot indicate cavities of varying 
size ri = 7io/o" + 0.5: 1.0 (circles), 1.5 (squares), 2.0 (left- 
triangles), 3.0 (right-triangles), 5.5 (up-triangles). Linear re- 
sponse approximation [Eqs. and (|14p ] was used to cal- 
culate Ec/Eq. The solid and dashed lines refer to continuum 
limits in Eqs. (|10|) and ([T}, respectively. 

creates surface charges. This picture is what we know as 
macroscopic dielectric described by material Maxwell's 
equations for which any field within polarized dielectric 
depends on the global shape of the sample This 
phase, which can be described as conventional dielectric, 
is shown by the shaded area in the space of parameters 
{e,i?i} in Fig. O For the parameters in the un-shaded 
area (marked as "polar liquid" in Fig. [2|) , the long-range 
polarization wave does not exist and any polarization 
wave in the liquid decays on a microscopic length Q. 
The local field is then independent of the sample shape 
and the rules of macroscopic electrostatics do not apply. 
This regime of relatively low polarities has an approxi- 
mate solution given by continuum limit of Eq. (jTU]). 

The quasi-continuum result of Eq. pop in fact corre- 
sponds to the separation of length-scales in which the 
cavity radius is larger than the correlation lengths for 
both the longitudinal and transverse dipolar correlations, 
i?i > Al,At. The use of Eq. (O then largely elimi- 
nates the effect of the transverse response on the con- 
tinuum portion of the response function. The transverse 
correlation length At is an increasing function of sol- 
vent dielectric constant (l3| growing to infinity at the 
ferroelectric transition. Therefore, the semi-continuum 
limit should become invalid at some e, and that happens 
through a discontinuous branching of the continuum so- 
lution between Eqs. pO)) and ([1]). In order for a solution 
to switch to the ordinary macroscopic limit, the singu- 
larity k* should be a part of the sample's spectrum of 
wavenumbers. The spectrum of k is limited to a dis- 
crete set of lattice values for a finite-size sample, and it 
is hardly possible for k* to coincide with one of the lat- 
tice vectors. Indeed, when continuous /c-integration in 
Eq. d?]) is replaced with the lattice sum according to the 
rule 

I dk/{2n f ^ J2 ' (12) 
we do not observe a rising part of the cavity field (dash- 



dotted fine in Fig. [T]). In Eq. (fT^ . L is the size of 
the cubic lattice and the lattice wavevectors are k = 
{2Tr/ L){n, I, m}. As expected from this calculation, we 
in fact have not observed switching to the ordinary con- 
tinuum in our numerical simulations. 

We have carried out MC simulations with the stan- 
dard NVT Metropolis algorithm, periodic boundary con- 
ditions, and the cutoff of the dipolar forces at L/2 (see 
Ref.lB for the details of the simulation protocol) . The 
initial configuration was set up as face-centered cubic lat- 
tice with random dipolar orientations and varied number 
of particles N. The spherical cavity was created at the 
center of the simulation box and the solvent diameter was 
adjusted to produce the bulk density p* — pa^ — 0.8. 
Reaction-field corrections with the dielectric constant 
equal to that of the liquid (from separate MC simula- 
tions) were used for the dipolar interactions to speed up 
the simulations. The results were identical to simulations 
employing Ewald sums. 

The cavity field was calculated from the linear response 
approximation according to the equation: 

EjEo - 1 + (/3/3)((5E, • 5M.) - E,,,,. (13) 

Here SEs and (SM are the fluctuations of the field at the 
cavity center and the total dipole moment of the liquid, 
respectively. The term Scom derived here from the pro- 
cedure suggested in Ref. [3, corrects for the cutoff of the 
dipolar interactions at the distance Tc in the simulation 
protocol: 

2(e-l) / £-1 fRiV\ , , 

Figure [3] shows the cavity field from MC simulations. 
The predictions of two continuum solutions, Eqs. H]) and 
(fTU)) . are shown by the dashed and solid lines, respec- 
tively. It turns out that both qualitative predictions of 
continuum electrostatics [Eq. ([1])] are violated. First, 
there is a significant dependence of E^ on the cavity size, 
as expected for cavities comparable in size to the liquid 
particles. Second, the simulated dependence Ec{e) never 
reaches the continuum limit of Eq. ([T]) , but instead levels 
off with increasing cavity size at the solution given by 

Eq. nni). 

The upward deviation of simulated cavity fields from 
Eq. pTI)l is a consequence of a particular orientational 
structure on the cavity's surface. Figure H] shows the 
distance dependence of the orientational order parameter 
formed by projecting the unit dipole vector, e^, on the 
unit radius vector, fj = Yj/rj: 

p{r)^(^P2{v,-e,)5{v,-v)^, (15) 

where P2{x) is the second Legendre polynomial. The 
surface dipoles tend to orient orthogonally to the surface 
normal with increasing polarity, a behavior well docu- 
mented for 2D dipolar fluids 16]. Surface orientation of 
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trostatic free energy of polarizing the dielectric [l| does 
not change at the branching point: 
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FIG. 4; Orientational order parameter p(r) [Eq. p5|l ] of the 
dipoles surrounding the cavity {Ri/a = 1.0) for varying e: 
1.4 (up-triangles) , 8.5 (down-triangles), 17 (diamonds), 30.6 
(squares), 53.7 (circles). The inset shows the splitting of the 
correlator {SEa ■ SNL) in Eq. p3|l into contributions from the 
first solvation shell (dotted line), from the second solvation 
shell (dashed line) , and first and second shell combined (dash- 
dotted line). The solid line indicates the overall correlator. 



dipoles leads to overscreening of the external field such 
that the electric field from the first solvation shell is di- 
rected oppositely to the external field (inset in Fig. [5]). 
This effect is partially compensated by a positive field 
from the second solvation shell, and it takes several shells 
to make the overall cavity field. For larger cavities (not 
shown here), the field of the first two solvation shells 
makes almost the entire cavity field such that the solvent 
response is more local and continuum-like. 

In conclusion, we have followed the procedure, first 
suggested by Maxwell [17], to measure microscopic fields 
in a polarized dielectric by carving cavities in it. The 
combination of numerical simulations and analytical the- 
ory drew a new picture of what is commonly called a 
macroscopic dielectric. We found that fields of macro- 
scopic electrostatics are formed only for sufficiently large 
polarity of the liquid and the cavity size as a singularity 
in the microscopic response function producing a non- 
decaying longitudinal polarization wave. The total elec- 
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where M is the total dipole of the polarized liquid. 
The integrated dipole moment of the cavity Mc — 
—3floP/{2e+l) does not depend on which solution for 
the cavity field is realized. Therefore, the decrease in 
the electrostatic energy of the cavity, caused by a lower 
cavity field, is released to the longitudinal wave. The 
appearance of this solution within the theory depends 
on the order of continuum {Rq/{A]^,At) ^ 1) and ther- 
modynamic (L — > oo in Eq. p2p ) limits. It is up to 
experimental measurements of cavity fields to determine 
which limit should be taken first. 

The cavity size reached in our simulations, 2Rq ~ 5 
nm, is of the order of that commonly realized for small 
nanoparticles, given the typical length-scale of molecular 
liquids cr ~ 4 — 5 A. We could never reach the limit of 
macroscopic electrostatics on that length-scale suggesting 
that electrostatics of nanocavities is not consistent with 
material Maxwell's equations. Macroscopic electrostatics 
[l[ assumes that polarization is a continuous field termi- 
nated at the interface where it creates a surface charge. 
The density of the surface charge is equal to the polar- 
ization projection normal to the surface 0. Restructur- 
ing of the liquid interface, in particular eliminating the 
normal polarization projection (Fig. |4|), will eventually 
affect the fields within dielectrics, cavity fields included. 
We need to notice that no direct measurements of fields 
within microscopic cavities in polar liquids have been, to 
our knowledge, reported in the literature. Experimen- 
tal evidence may arrive from measurements of dielectric 
relaxation of photoexcited dipolar impurities. 
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